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ABSTRACT 

Galaxies inhabiting a cluster environment experience significant evolution in their 
orbital motions throughout time; this is accompanied by changes in the anisotropy 
parameter, measuring the relative importance of radial and tangential motions for a 
given class of objects. Along with orbital changes, galaxies in clusters are well known to 
undergo severe alteration in their hot/cold gas content and star formation properties. 
Understanding the link between the changes in the internal properties of galaxies and 
their orbital motion is of crucial importance in the study of galaxy evolution, as it 
could unveil the primary mechanism responsible for its environmental dependence. 
Do the changes in the internal properties happen in parallel with those in the orbital 
motion? Or are the orbital features at the time of infall what determines the fate 
of the member galaxies? Alternatively: are the properties of galaxies at a given time 
related to the coeval orbital anisotropy or are they better related to the anisotropy at 
infall? In order to answer these questions, we studied the orbital evolution of different 
galaxy populations in the semi-analytic models of Guo ct al. (2011) applied on to the 
Millennium Simulation. For each class of objects, characterised by different internal 
properties (such as age, star formation rate and colour), we studied the anisotropy 
profile at redshift zero and its evolution by tracing the progenitors back in time. 
We conclude that the orbital properties at infall strongly influence the subsequent 
evolution of the internal features of galaxies and that the overall anisotropy of the 
galaxy population tends to increase with time. 

Key words: methods: numerical - galaxies: evolution - galaxies: kinematics and 
dynamics - cosmology: theory, large-scale structure of Universe. 



1 INTRODUCTION 

In the currently favoured scenario for the formation of cos- 
mological structure in the Universe, dark matter halos form 
and merge giving rise to a hierarchy of objects assem- 
bling in a bottom- up fashion. In this context, galaxy for- 
mation is pictured to arise from cooling and condensation of 
baryons within the potential well associ ated to these dark- 
matter structures (|White fc Reeslll978l ). Once the galaxy 
is formed, its evolution will be driven by (i) "nature", 
i.e. the object's intrinsic features (essentially stellar mass) 
and (ii) "nurture", external processes related to the envi- 
ronment the galaxy inhabits during different stages of its 
history. It is indeed well known that, although the struc- 
tural properti es of galaxies are mainly determined by their 
stellar mass ( Kauffma nn et alj 120031 ; iTanaka et al. I |2004 
Ivan den Bosch et al.l 20081 ). the existence of an environmen- 
tal dependence ca n not be disregarded llHogg et al. 20041: 
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120051 , among the others) . A combination of these effects is re- 
sponsi ble for the observed "bimodality" in the gal axy prop- 
erties (|Baldrv et al.l 12004 IKauffmann et al.l 120041 ). namely 
the existence of two well-distinguished classes of objects 
characterised by either red colour/high mass/old stellar pop- 
ulation or blue colour/lower masses/young stellar popula- 
tion, th e former resid i ng preferential l y in over-dense environ- 



ments dOemledll974lDrcsslcr 198p|; iBower fc Baloghl 



Balogh et alJ |2004 iBall et al.1 120081 : iBamford et all 



2004; 



Skibbaet al.Tl2009r i. while the latter being found mainly 



2009; 



the field. A clear and thorough physical picture of galaxy 
evolution is yet to come; understanding what mechanisms 
play the leading role in shaping galaxy properties is a 
topic of ever increasing interest and research activity, but 
the results are still controversial. Among the environmental 
processes, four broad classes of mechanisms are generally 
considered: galaxy mergers, negligible in massive clusters, 
but important in small groups, which drive morphological 
changes and can affect star formation (IToomre fc Toomre 



changes a nd can attect star to rmation (lloomre fe lo omr( 
19721 : iFarouki fc Shapirdll98ll ): strangulation ( Larson et al 



1980D . a sum of effects leading to the removal of the hot 
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gas halo associated to a galaxy when this is accreted onto 
a lar ger structure; ram-pressure stripping (|Gunn fc Got j 
1972), important in massive clusters where the density of 
the intra-cluster medium is highest, it leads to the progres- 
sive stripping of the gaseous component (both hot and cold) 
of the satellite galaxies; tidal effects, arising from the gravi- 
tational interaction with other members and with the cluster 
potential itself, they cau se stripping and heating jRichstond 
1 19761 ; iMoore et alii 1996ft . Some of these processes are bet- 
ter understood and formalised (e.g. ram-pressure stripping), 
while others still lack a strict, physical description and are 
referred to in more generic terms (e.g. strangulation); on 
top of this, the regimes where each of the mechanisms is 
active/ unimportant are only broadly assessed. In galaxy 
clusters, the densest possible environments, the contribu- 
tion of mergers can be neglected, due to the large veloc- 
ity dispersion of the system; strangulation occurs rapidly 
and devoids the hot-gas reservoir, while ram-pressure strip- 
ping and tidal interactions proceed as the satellite plunges 
into its host. The importance of these last two processes 
increases with local density (gas density for ram-pressure 
stripping, total matter density for tidal effects), but while 
in the first case this results into enhanced stripping, the ef- 
fect of tidal interactions is believed to be mainly that of 
an induced ga s consumption following the increase of nu- 
clear activity (|Boselli fc Gavazzj 12006,1. Given the depen- 
dence upon gas and matter density and how these are, in 
turn, related to the radial distance from the cluster centre, 
it is obvious to expect these processes to strongly depend 
on the orbital history of the satellites, namely how often 
these happen to transit the innermost regions with respect 
to the outskirts. This will depend on the initial orbit of the 
galaxy at the time of accretion, as well as on its time evo- 
lution within the dynamically-active cluster region. There 
exist observational evidences that Hi-deficient galaxies in 
nearby clusters are early-type spirals on radial orbits, while 
gas-ri ch galaxies are characterised by m ore tangential mo- 
tions (|Dresslerjll986l ; ISolanes et al.ll200lh ; other works stud- 
ied the differences in the velocity distribution of early-type 
and late-type galaxies and found evidences of the latter mov- 
ing on slightly more radial orbits, especially at large clus- 
tercentric radii dMahdavi et al.l Il999l; Ivan der Marel et al 
200d; iKatgert et al.ll2004l; iBiviano fe Katgertll2004l; iBiviano 



2009 , 120081 ; IBiviano fc Poggiantil 120091 ; IWoitak fc Lokas 
20101 '). Deriving this information from a sample of observed 
cluster members is not an easy task and various techniques 
were developed appositely. By instance, knowing the pro- 
jected number density profile and line-of-sight velocity dis- 
persion of the selec ted galaxies, one could pe rform a Jeans 
dynamical analysis |Binnev fc Tremaine!ll987M to derive the 
cluster mass and velocity anisotropy profiles (quantifying 
the relative importance of radial and tangential orbits as 
a function of the clustercentric radius). This technique re- 
lies on the assumption of spherical symmetry, collisionless 
dynamics and dynamical equilibrium, in addition to be- 
ing hampered by the so-called "mass-anisotropy" degener- 
acy (even if co nsiderable progress has recently been done to 



remove it, see 



Lokas fc Mamonll2003l; iBattaglia et aI.ll200Sl ; 



evolution of galaxies is linked to their specific orbital mo- 
tion. Is the type of orbit determining the efficiency of the 
environmental processes in shaping the internal properties 
of the object? In this case we would expect severely affected 
galaxies to move on more markedly radial orbits, as these 
cross the cluster right down to its innermost regions. Or do 
the changes in the galaxy properties happen in parallel to 
those in the orbital motion? In this case we would expect the 
orbits of the objects having suffered major environmental in- 
fluences to differentiate from those of just-accreted satellites. 
The importance of the numerical approach in the study of 
galaxy evolution need not be stressed; cosmological sim- 
ulations provide a fundamental tool to understanding the 
assembly of structures throughout time, but the reliability 



ity is subject of debate 


(see, e.g., ScannaDieco et al. 


12011). 


Semi-analytic models ( 


White & Frenkl 199ll; Cole 


1991; 


Kauffmann et al.l Il993l; 


Cole et al. 19941; Kauffmann et al.l 


1999]; Suringel et al. 200ll; see Baughl 20061 for a review) 



IWoitak et al.l2009l ). This limitations, along with the difficul 
ties in performing observations in the outskirts of clusters, 
should induce caution in the interpretation of the results. 
Yet it would be interesting to know the extent to which the 



present a powerful, hybrid approach to the problem of galaxy 
formation and evolution: they make use of accurate, dark- 
matter-only simulations to account for the growth of struc- 
tures in the cosmological context and regulate the formation 
and evolution of galaxies according to a set of analytical pre- 
scriptions encompassing the physics of the baryonic compo- 
nent. In this work we use one of the most advanced semi- 
analytic models developed to date to study the link between 
orbital and internal properties of galaxies belonging to mas- 
sive clusters in a ACDM cosmology. We will show how the 
results of this analysis suggest a scenario where the specific 
orbital features of the satellites at the time of infall have ma- 
jor consequences on their evolution. The paper is organised 
as follows: Sec.[5]describes the simulation, the semi-analytic 
model and the choice of the cluster sample used in our analy- 
sis; Sec.[3]reports the results for the anisotropy parameter of 
the member galaxies at redshift zero (|3.1|l . higher redshifts 
(|3.2|l , at the time of the last infall inside their host (|3.3[) and 
regarding its time evolution (|3.4p ; finally, Sections 2] and [S] 
contain a brief summary and discussion of the results. 



2 THE SEMI-ANALYTIC MODELS AND THE 
SELECTED SAMPLE 

The samples analysed in this work were extract ed from the 
galaxy catalogue obtained by iGuo et all (|201lh (hereafter 
GUOll) after running their galaxy for mation models on th e 
Millennium Simulation (hereafter MS; ISpringel et al.ll^OoBT ). 
The MS shows how dark matter structures form and evolve 
in a ACDM cosmological scenario characterised by the pa- 
rameters fitot =1, fim =0.25, Q, b = 0.045, fi A = 0.75, h = 
0.73, <T8 = 0.9, n s — 1. The evolution of cosmological struc- 
tures is traced by 2160 3 particles moving in a periodic box 
of side 500 Mpc/h under the mutual gravitational influence. 
The results of the simulation were stored at 64 different 
times starting from redshift 127 down to redshift zero. At 
each of these times a catalogue of bound structures and 
substructu res was generated , by applying a friend-of-friend 
technique jDavis et al1ll985l ) and the SUBFIND algorithm 
(|Springel et al.ll200ll ) on the particle distribution; these cat- 
alogues constitute the basis for recovering the merger history 
of structures throughout time, also referred to as the merger 
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tree of the simulation. Semi- analytic models of galaxy for- 
mation, as those of GUOll, provide a description of the 
cosmological evolution of baryonic matter; a gas distribu- 
tion is associated to each bound structure identified in the 
simulation and its evolution in time is regulated by a set of 
recipes implemented on to the dark-matter merger tree iden- 
tified in the simulation: according to the specific histories of 
each of the substructures, the evolution of the associated 
baryonic components will follow its own, peculiar path. The 
physical processes implemented in the recipes of GUOll in- 
clude cooling, star formation, supernova and AGN feedback, 
hot-gas stripping, metal enrichment and alone they provide 
a remarkable match to the abundance and large-scale clus- 
tering properties of the observed galaxy population at low 
redshift. We will comment on the limits of the model in the 
discussion of our results; we refer the reader to GUOll for a 
thorough description of the implemented physics and of the 
strengths and weaknesses of their approach. 
We have selected the 1000 most massive clusters identified 
in the MS at redshift zero; these have virial massed greater 
than 2 x 10 14 M@ and contain a total of around one mil- 
lion galaxies. Although all of these galaxies are originally 
associated to a specific dark-matter subhalo, not all of them 
preserve their dark-matter component throughout the sim- 
ulation time; this means, in turns, that in these cases the 
dynamical evolution can no longer be provided by the under- 
lying simulation. The fate of these orphan galaxies (gener- 
ally referred to as "type 2s" in the galaxy formation model), 
which have lost their dark-matter component, is to eventu- 
ally merge with the galaxy at the centre of the cluster they 
inhabit, within the timescales set out by dynamical friction. 
During this time, the orbit of the orphan galaxy is traced by 
the most bound particle present in its dark matter halo be- 
fore this vanished, modified by a shrinking factor introduced 
to mimic the orbital decay caused by dynamical friction. 
Since the orbits of these galaxies are altered, we decided not 
to include them in our analysis. This leaves us with « 2 x 10 5 
galaxies in the mass range 10 3 < M*/M Q < 10 12 , where M„ 
indicates stellar mass. We have stacked the resulting galaxy 
sample by subtracting bulk motions and by normalising po- 
sitions and velocities to the virial values. To summarise, we 
are studying the mean behaviour of galaxies residing in the 
most massive clusters at redshift zero. 

From the galaxy catalogue we have access to a wealth of 
information regarding the internal properties of the selected 
objects; those we are primarily interested in are stellar/dark 
matter masses and colour, the latter defined as the the differ- 
ence between the rest-frame total absolute magnitude in the 
SDSS u and i bands. We will eventually discuss the impact 
on our results of using other quantities, such as mean stel- 
lar age and specific star formation rate, to split our galaxy 
sample into distinct populations. 



3 RESULTS 

In this section we will present the results on the velocity 
anisotropy for the galaxies in the selected sample. This quan- 
tity is defined as: 



= 1- 



(1) 



where a t and a r are the velocity dispersions in the tangential 
and radial directiorQ respectively. The velocity anisotropy 
can take values ranging from -co to 1; the former case cor- 
responds to purely circular motions (no dispersion in the ra- 
dial direction), while in the latter the orbits are completely 
stretched along the radial direction (no tangential disper- 
sion) . The case /? « is referred to as isotropy and it occurs 
when the velocity dispersion is of comparable magnitude in 
both radial and tangential directions. We will mostly dis- 
cuss global values of /3, describing a full population of ob- 
jects without discriminating the spatial distribution of the 
members; in one case we will also show the radial profile, to 
underline both the general properties of the anisotropy pa- 
rameter and its dependence on the distance form the centre 
of the cluster. In what follows, we split the galaxy popula- 
tion into a red and blue sample according to the value of the 
u — i colour indicator, namely whether it is greater or smaller 
than a certain threshold; this is set to the value at which the 
colour distribution of the galaxy sample splits into two well- 
defined, different components. At redshift zero, this happens 
around u — i = 2.5; considering all the galaxies within three 
virial radii, the resulting "red" sample consists of « 72000 
objects, against the fs 130000 of the "blue" counterpart- 



3.1 Anisotropy at redshift zero 

We first analyse the results at redshift zero and as a func- 
tion of the radial distance of the galaxies from the centre 
of the stacked cluster. The curves in Fig. [T] show the be- 
haviour of the anisotropy parameter in 20 equally-spaced 
radial shells, extending from roughly 0.1 to 3 virial radii. 
The error bars were calculated by means of a bootstrapping 
algorithm and correspond to two standard deviations (this 
holds for all the uncertainties reported in this paper, un- 
less explicitly stated otherwise). The black curve shows the 
result for the full population: the importance of radial mo- 
tions increases, moving from the central regions outwards 
and peaks around 1.5 — 2 virial radii. This beha viour is typi- 
cal o f systems formed by gravitational collapse (|van Albadal 
Il982h and it i s compatible with othe r results in the literature 



Sales et al. 2007 


; Woitak et al. 20091: Biviano & Poeeianti 


20091: Host et al. 


120091: Lemze et al.l 2011: LaDi & Cavaliers 



201 if ) . The blue and red curves represent the results for the 
subgroups of galaxies characterised by blue and red colours, 
respectively; the clear message coming from the plot is that 



These are defined as the masses of the spherical regions centred 
on the potential minimum of the smooth halo and corresponding 
to an overdensity of A, typically as 200, with respect to either the 
critical density p cr ;t or to the mean background density p m = 
Qmpcrit- Hereafter, when referring to either the virial radius or 
the virial mass of an object, we assume the overdensity to be 
defined with respect to the critical density. 



2 By tangential direction we mean that identified by variations 
in any of the angles 9 or <j>, i.e. <r 2 = (cr 2 + cr 2 )/2. 

3 The higher fraction of blue objects can be traced back to (i) the 
exclusion of type 2s and (ii) the choice of three virial radii. When 
restricting the analysis to 1.5 virial radii and including type 2s, 
the number of red objects becomes nearly four times as large as 
that of the blue ones. 
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Anisotropy profile 
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Figure 1. Anisotropy profile for the stacked sample at rcdshift 
zero. The black, solid curve shows the radial behaviour of the 
anisotropy parameter for the full population of galaxies found in 
the selected cluster sample; the red (blue), solid curve represents 
the result for the subsample of galaxies characterised by an u — i 
colour greater (smaller) than 2.5. The error bars are evaluated 
via a bootstrapping algorithm and correspond to two standard 
deviations. The horizontal, dotted lines mark the global value of 
the anisotropy parameter for all the galaxies within three virial 
radii. 

the blue population has a systematically and significantly 
lower anisotropy than the red population, at least between 
0.5 and 2.5 virial radii. This is confirmed by the global val- 
ues of /3, marked by the horizontal, dotted lines. These were 
computed out of all the galaxies found at distances less 
than three virial radii from the centre, with the following 
result: p aU = 0.253 ± 0.006 ,/3 Wtle = 0.189 ± 0.006 , /3 red = 
0.345 ± 0.008. In summary, red galaxies in our sample move 
on more radially-biased orbits than the blue counterparts. 
A clear explanation of this effect will progressively arise in 
the following sections. 

3.2 Anisotropy at high redshift 

In what follows, we focus our attention on the time evolu- 
tion of the anisotropy parameter. Having not found signif- 
icant changes in the radial behaviour with respect to the 
trend shown in Fig. [T] we will only refer to global values of 
/3 hereafter. We proceeded in two ways: we considered all 
the member galaxies belonging to the high-redshift progeni- 
tors of the selected clusters and, in parallel, we also analysed 
the progenitors of the redshift-zero galaxies only. The two 
approaches differ in that the first galaxy sample may con- 
tain objects that do not survive until redshift zero, whereas 
in the second case only satellites that have a redshift-zero 
descendant are considered. We show the results of the first 
method in Fig. [21 where the value of /3, computed for the 
full population of member galaxies regardless of their spatial 
position, is plotted as a function of redshift. The colour code 
is the same as for Fig. [TJ keeping the threshold u — i — 2.5 
for the splitting between the red and blue population or ad- 
justing it to the variation in the colour distribution of the 
high-redshift galaxies introduces no substantial changes in 
the results. 

The two main considerations arising from this analysis are 
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Figure 2. Global value of the anisotropy parameter as a function 
of redshift. Considered are all the member galaxies belonging to 
the high-redshift progenitors of the redshift-zero cluster sample; 
these may include objects which do not survive to redshift zero. 
The black curve shows the result for the full galaxy population at 
each redshift, whereas the red (blue) curve refer to the subsample 
of objects with u — i colours greater (smaller) than a threshold; 
this is set to be the value at which the colour distribution splits 
into two, well-defined component (e.g. u — i PS 2.5 at redshift 
zero, u — i « 2.35 at z = 0.7). The error bars are evaluated 
via a bootstrapping algorithm and correspond to two standard 
deviations. 

(i) that the global value of the anisotropy parameter does 
not significantly evolve in time and (it) that the blue pop- 
ulation is characterised by a lower degree of anisotropy at 
all redshifts. A more interesting result is given by the analy- 
sis of the second sample of high-redshift galaxies. These are 
progenitors of redshift-zero objects that are already mem- 
bers of their final, host cluster at the redshift of interest 
and that will not leave it anytime afterwards. Fig. [3] shows 
the global anisotropy for this class of galaxies; the colour 
code refers to the redshift-zero population only and not to 
the progenitors: the red (blue) curve corresponds to progen- 
itors of galaxies that are red (blue) at redshift zero. The 
plot shows with striking clarity that, at each redshift, the 
anisotropy of the progenitors of redshift-zero galaxies is con- 
siderably lower than the value for the whole, high-redshift 
population (Fig. [2] overplotted in gray); this effect seems to 
be even stronger for progenitors of galaxies that are blue at 
redshift zero. These results suggest (i) that galaxies on more 
radial orbits have less chances to survive and are progres- 
sively removed from the member population and (ii) that 
a similar selection effect acts on the galaxy colour (namely: 
galaxies on radial orbits have very low chances to remain 
blue until redshift zero). The second point will become even 
clearer when analysing the anisotropy at the time the galax- 
ies fall inside their host cluster. 



3.3 Anisotropy at infall 

We finally study the orbital properties of galaxies at the 
time they become members of their cluster; by this we mean 
the moment the dark-matter halo of a galaxy is not an in- 
dependent structure anymore, but a substructure identified 
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Figure 3. Global value of the anisotropy parameter as a function 
of rcdshift. Considered are the progenitors of the galaxies belong- 
ing to the redshift-zero cluster sample; only progenitors which are 
already satellites and will remain satellites down to rcdshift zero 
arc taken into consideration. The black curve shows the result for 
the full population at each redshift, whereas the red (blue) curve 
refer to the subsamplc of objects with redshift-zero descendant 
characterised by u — i colours greater (smaller) than 2.5. Over- 
plotted in gray are the results displayed in Fig. [2] The error bars 
are evaluated via a bootstrapping algorithm and correspond to 
two standard deviations. 

within a larger halo. We want to see how the anisotropy of 
the infalling population varies with time. Again, we split the 
analysis into two parts; first, we consider the full sample of 
galaxies joining the high-redshift progenitors of our original 
cluster sample and, second, we also restrict the study to the 
subsample of galaxies with a descendant at redshift zero. 
The results of the first approach are displayed in Fig. [4j 
where we plot (3 against the infall redshift. The red (blue) 
curve refers again to galaxies that are red (blue) at the red- 
shift under consideration; keeping u — i — 2.5, or varying the 
u — i colour threshold to adapt to the colour evolution of the 
galaxy population at high-redshift introduces non substan- 
tial changes in the results. A clear trend appears, showing 
that the infall anisotropy increases going towards lower red- 
shifts; this is not surprising, as mass-accretion is expected 
to occur progressively along small filaments, extending ra- 
dially outside massive clusters. The second feature which is 
apparent in the plot is that most of the infalling population 
consists of blue galaxies, as shown by the vicinity of the black 
and blue curve, as well as by the size of the error bars. Third 
we see that, again, the anisotropy of blue galaxies is lower 
than that of the red galaxies. At a first glance, this may re- 
sult somewhat unexpected as there is no obvious reason why 
such a trend should already be in place before entering the 
cluster environment. We note, though, that the mean mass 
of the red sample is considerably higher than for the blue 
counterpart (almost two order of magnitudes for both stel- 
lar and dark-matter masses); previous works have already 
stated that more massive satellites tend to approach the 
cluster along more eccentric orbits, as they are more likely 
to re side in dense, radial filaments than the least massive 
ones (|Tormenll 19971 ). This effect could explain the behaviour 
observed in Fig. [4] and also play a role in the interpretation 
of Fig. [T] an initially higher anisotropy for the red popula- 
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Figure 4. Global value of the anisotropy parameter for the in- 
falling population, as a function of infall rcdshift. Considered are 
all the galaxies falling inside the progenitors of the redshift-zero 
clusters at different times; these may include objects which do not 
survive to redshift zero. The black curve shows the result for the 
full infalling population at each redshift, whereas the red (blue) 
curve refer to the subsample of objects with u — i colours greater 
(smaller) than 2.5. The error bars are evaluated via a bootstrap- 
ping algorithm and correspond to two standard deviations. 

tion of infalling galaxies may leave an imprint in the final 
profile at redshift zero. The impact of this effect is anyway 
limited by the small number of these galaxies within the in- 
falling population: at z = 0.7, out of 19049 infalling galaxies 
only 419 are red. Clearly, these objects may only represent 
a contribution to the population of red, member galaxies at 
low redshifts. Fig. [5] shows the results for the subsample of 
infalling galaxies with a redshift-zero descendant. As usual, 
in black are the result for the full population whereas the 
red (blue) curve refers to the subsample of galaxies that are 
progenitors of objects characterised by a red (blue) colour 
at redshift zero. Again we see that, at each redshift, the 
anisotropy of these galaxies is systematically lower than the 
value for the full, infalling population (Fig. [4] overplotted 
in gray). This suggests that galaxies present in the clus- 
ter at redshift zero tend to originate from the subsample of 
objects entering the cluster environment with the least ra- 
dially stretched orbit. We then split the full population not 
according to the colour of the redshift-zero descendant, but 
on the basis of their colour at the time of infall; the orange 
(cyan) curve represent the infall anisotropy of galaxies that 
are red (blue) at the time of infall. We clearly see that the 
galaxies which are blue at redshift zero are descendant of 
the subgroup of infalling, blue galaxies with the lowest in- 
fall anisotropy. It seems, therefore, that not only does the 
initial orbit strongly influence how long the galaxy is going 
to survive in a cluster environment, but that it also plays an 
important role in determining the evolution of its internal 
properties. 

3.4 Evolution of anisotropy in time 

When visually comparing the results for the infall anisotropy 
(Fig. [4] and [5} to the typical values at redshift zero (Fig. [T] 
and first points of Fig. [2} we already guess that p increases, 
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Figure 5. Global value of the anisotropy parameter for the in- 
falling population, as a function of infall rcdshift. Considered are 
the progenitors of the galaxies belonging to the redshift-zero clus- 
ter sample; only progenitors which are becoming satellites for 
the last time in their history are taken into consideration. The 
black curve shows the result for the full infalling population at 
each rcdshift, whereas the red (blue) curve refer to the subsamplc 
of objects with redshift-zero descendant characterised by u — i 
colours greater (smaller) than 2.5. Displayed in orange (cyan) is 
the anisotropy of the subgroups of infalling satellites which are 
red (blue) at the time of infall. Ovcrplottcd in gray are the results 
displayed in Fig. [4] The error bars are evaluated via a bootstrap- 
ping algorithm and correspond to two standard deviations. 

from the time the galaxies enter the cluster to the end of 
the simulation. This is shown more quantitatively in Fig. [6] 
Each of the curves represents the results at a specific red- 
shift, from 0.7 (in red) to zero (in black); the points on each 
curve correspond to the anisotropy, at this redshift, of the 
subgroup of galaxies that fell inside the clusters at the red- 
shift given in abscissa. The first points of the lines show the 
anisotropy of the galaxies which are just infalling and, taken 
altogether, they reproduce the black curve of Fig. [4] Overall, 
there is a strong indication of an increase of j3 with time, 
after the galaxies become satellites and orbit within the en- 
vironment of a larger halo; this seems to occur rapidly within 
the first 2 Gyrs after infall and more gradually afterwards. 
In order to understand what this increase was due to, we 
investigated the evolution of individual orbits. From a ran- 
dom selection of 500 galaxies belonging to the redshift-zero 
sample, we isolated ~ 100 objects sharing the same in- 
fall snapshot and have recovered their full orbital history. 
Knowing the satellites' mass, their position, velocity and the 
virial mass of the hosts at the time of infall, we integrated 
the orbits forward in time. We used the leapfrog method 
to solve for the motion of the objects, with timesteps of 
around 90 Myrs. This time resolution corresponds, approx- 
imately, to one fourth of the temporal separation between 
the different snapshots in the MS (reducing the size of the 
timestep further does not significantly change the results). 
The integration proceeds for a number of steps correspond- 
ing to the total time elapsing from infall [z — 1) to red- 
shift zero. We adopt two different approaches in the integra- 
tion: in the first case (hereafter case A) we keep the mass of 
both the host and the satellite constant at the infall value, 
whereas in the second case (hereafter case B) we update 



both masses according to the values provided in each of the 
snapshot^] and by means of a linear interpolation to ac- 
count for our finer temporal integration. In both cases we 
assu me the mass of the host to follo w a spherical NFW pro- 
file (|Navarro. Frenk.~fc Whitejll996l ). The concentration pa- 
rameters were evaluated from Eq. 5 of iNeto et alj (|2007h . 
who obtained the concentration-mass relation for the halos 
identified in the MS at redshift zero; we do also allow the 
concentration parameter to vary with redsh ift and assume 
the dependency found by iDuffv et al.l (|200ct ) for simulated 
halos in the redshift range — 2. In both integrations the 
satellites are regarded as point masses and no treatment of 
dynamical friction is included; this sh ould not significantly 
impa c t the shape of the orbits, tho ugh (|van den Bosch et alj 
1 19991 ; iPenarrubia fc Benson! [20051 ). At each timestep, the 
mass felt by the orbiting satellites varies according to the 
distance from the centre of the host (for both case A and 
B) and to the cosmological accretion (taken into account in 
case B only) . Case A hardly ever reproduces an orbit close to 
the original, whereas case B provides a better description to 
the motions occurring in the simulation. Even in this case, 
however, a perfect match to the original orbit is not guar- 
anteed. Our approximation for the distribution of the host 
mass, on top of the fact that we are neglecting the effect 
of interactions with other orbiting galaxies, can sometimes 
result in evident mismatched]. Fig. [7] shows few examples of 
original orbits from the MS (black curves in the plots on the 
left column) and the results of our integration for case A and 
B (overplotted purple and magenta curves, respectively); as 
mentioned, in some cases our approximation grasps the orig- 
inal dynamics pretty well (e.g. magenta curves in the first 
three rows) , while in others we are clearly missing something 
(e.g. last row). The right column shows the time evolution 
of the mass felt by the satellite in both case A and B and 
explains the differences in the corresponding orbital evolu- 
tion. 

Notwithstanding the simplicity of our model, we can repro- 
duce fairly well the evolution of the anisotropy parameter 
for the selected galaxies. This is shown in Fig. [5] where j3 
is plotted against the snapshot number, the latter ranging 
from 41 (i.e. z = 1, the infall redshift) to 63 (i.e. z — 0, 
the end of the simulation). The green curve shows the result 
obtained from the simulation catalogue, i.e. the real tem- 
poral evolution of the anisotropy parameter for the selected 
galaxy sample. In purple (magenta) are the results from case 
A (B). We have split the integration into five slots lasting 
each around 1 Gyrs; at the end of each slot, we updated 
the position and velocity of the satellites from the real value 
obtained from the MS. In case A, /3 immediately increases 
towards one after the integration is switched on; this is due 
to most of the galaxies becoming unbound and can be re- 
lated to a progressive underestimation of the depth of the 

4 The value for the dark-matter mass of the satellite, recovered 
from the simulation, is subject to inaccuracies whos e importance 
depen ds on the spatial position of the object (see iKnebe et alj 
l201ll) . We have performed the integration updating the mass of 
the host only and keeping the mass of the satellite fixed at the 
infall value: the differences with the results obtained in case B are 
negligible. 

5 For a thorough a ccount of the difficul t ies in reproducing satel- 
lite orbits, see, e.g., iLux. Read, fc Lakj i20ld ). 
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potential well they are moving in. Case B, instead, provides 
a remarkable match to the original profile. The agreement 
is maintained to a reasonably good level also when the in- 
tegration is let proceed uninterrupted from infall to redshift 
zero. 

Interestingly, we found that this trend does not translate 
naively into the corresponding behaviour of the circularity 
distribution. Knowing position, velocity and masses one can 
straightforwardly compute the circularity parameter: 



Evolution of anisotropy in time 



L_ 

17 



(2) 




where L is the angular momentum of the orbit and L c is 
that of the circular orbit characterised by the same orbital 
energy E. The circularity parameter has finite values only 
for bound orbits, i.e. those with E < 0; the limit r\ = cor- 
responds to purely radial orbits, whereas 77 = 1 to perfectly 
circular motions. At each timestep of the integration, we 
computed the circularity parameter for each of the bound 
orbits and looked at its evolution with time. Fig. [5] shows 
the temporal evolution of the mean of the circularity dis- 
tribution, for both integration A and B. As in Fig. [S] the 
results are shown for five integration slots, at the beginning 
of which the position and velocity of the satellite is updated 
to the simulation value. The purple curves, corresponding to 
case A, are flat: the mean circularity does not evolve with 
time. The same holds for the magenta curves, corresponding 
to case B, which also show no net evolution, other than for 
a mild increase in the first slot. The slope is what is sig- 
nificant here, while the normalisation of the curves and the 
differences among them are not strictly meaningful. Indeed, 
the first point of each slot is computed out of position and 
velocities extracted from the simulation, while the potential 
is provided by our simplified model; this makes the overall 
evaluation of the circularities not fully consistent. When let- 
ting the integration run uninterrupted from infall to redshift 
zero, the results are again consistent with no net evolution 
of the mean circularity of the galaxy population for both in- 
tegrations. This analysis is performed on the bound subset 
of the galaxy orbits, while the results shown in Fig. [H] are 
obtained out of all the galaxies under consideration. When 
computing the evolution of the anisotropy parameter for the 
bound subsample, the results for case B are practically iden- 
tical (most of the orbits are bound, in this case); for case A, 
instead, the number of bound orbits is too small to derive 
significant prediction for the temporal evolution of /3. This 
little experiment seems to suggest that taking into account 
the evolution of the potential well of the host cluster, due 
to cosmological accretion, is enough to reproduce the trend 
in the anisotropy parameter as obtained from the MS for 
the sample of galaxies under consideration. This is possible 
notwithstanding the over-simplistic assumptions of a per- 
fectly spherically-symmetric, smooth distribution and of an 
isolated, point-mass satellite. At the same time, we do not 
register a corresponding decrease in the mean circularity of 
the bound orbits, as one would naively expect, but rather 
an absence of significant evolution. 



Figure 6. Global value of the anisotropy parameter at different 
times and as a function of infall redshift. Each curve represents the 
results at a specific redshift, as stated in the legend. The points 
on the curves correspond to the anisotropy of galaxies which have 
last become satellites at the redshift given on the x axis. The error 
bars are evaluated via a bootstrapping algorithm and correspond 
to two standard deviations. 

4 SUMMARY 

We have studied the evolution of the anisotropy parame- 
ter for galaxies orbiting within the most massive clusters 
extracted from the MS; we have related the value of this 
parameter to the internal properties of the galaxies, as pre- 
dicted by the semi-analytic models of GUOll. Our findings 
can be summarised as follows: 

• At redshift zero, blue galaxies move on less radial or- 
bits than red galaxies do (Pbiue ~ 0.19 vs. /3 rec j sa 0.35; see 
Fig.P). 

• At higher redshifts, progenitors of redshift-zero objects 
move on less radial orbits with respect to the full population 
of member galaxies (see Fig. [3J); this is particularly true for 
progenitors of galaxies that are blue at redshift zero. 

• The orbits of infalling galaxies become increasingly ra- 
dial going towards redshift zero (see Fig. 2| . 

• Of all the galaxies entering the cluster at a certain time, 
those that will survive until redshift zero are the subgroup 
characterised by the most tangential orbits at infall (see 
Fig.©. 

• Of all the blue galaxies having a redshift-zero descen- 
dant (either red or blue) and entering the cluster at a cer- 
tain time, those that remain blue until redshift zero are the 
subgroup of infalling objects on the most tangential orbits 
(compare the blue and cyan curves in Fig. 0. 

• The orbit of satellite galaxies become increasingly ra- 
dial in time, especially just after the infall inside the cluster 
environment (see Fig. [6}. 

• This increase can be reproduced by a simple integration 
of the galaxy orbit, taking into account the cosmological 
growth of the halo mass (see Sec. 13.41) . 

These findings and their statistical significance strongly 
suggest that the orbital features at infall have a major 
influence on the subsequent evolution of the galaxies, i.e. 
on their survival time within the cluster and on their turn- 
ing from being star-forming/blue to passive/red objects. 
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Evolution of the anlsotropy parameter: MS, case A and 



ntegration A vs 









Figure 7. Examples of integrated orbits. The left column shows 
the motion of four different galaxies around their host. The blue 
crosses represent the centre of the host and the initial position of 
the galaxy, while the circle gives the size of the virial region at 
the beginning of the integration. The black curve corresponds to 
the original orbit from the MS and the purple (magenta) curve 
shows the result of the integration in case A (B). Whenever the 
orbit is bound, a triangle is overplotted; this is not done for the 
original orbit, as the information on its energy is not available. 
The timestep is 90 Myrs, approximately one fourth of the tem- 
poral separation between two MS snapshots (sd 350 Myrs). The 
right column shows the evolution of the mass felt by the galaxy 
as it orbits around the host in both case A (purple curve) and 
case B (magenta curve); as spherical symmetry is assumed in the 
mass distribution of the host, this quantity corresponds to the 
mass contained within the sphere of radius r = r aa t — T*host- 



Galaxies on more markedly radial orbits will be more easily 
disrupted and removed from the cluster members as time 
goes by; even if they survive until redshift zero, these 
objects are those most likely to undergo changes in their 
internal properties and turn from being star-forming/blue 
to passive/red. Galaxies which, conversely, join their host 
with a significant component of their motion along the 
tangential direction tend to be less affected by the cluster 
environment; to this category belong all the blue galaxies 




50 55 
Snapshot 

Figure 8. Evolution of the anisotropy parameter as obtained 
from the MS (green curve) and for cases A (purple curve) and B 
(magenta curve). The integration of the orbits has been performed 
on a randomly selected subsample of approximately 100 galaxies 
sharing the same infall snapshot, given their initial position and 
velocity. The results of the integration are compared to the profile 
computed out of simulated data for the same galaxies. The plot 
shows the temporal evolution of /3 for the whole population, as 
a function of snapshot number - from 41 (i.e. z = 1, the infall 
redshift) to 63 (i.e. z = 0). The integration has been split into five 
slots extending roughly 1 Gyr in time; at the beginning of each 
slot (marked by a vertical, dotted line) the position and velocity 
of the satellite was updated to the real value obtained from the 
MS. The error bars are evaluated via a bootstrapping algorithm 
and correspond to one standard deviation. 



Evolution of the circularity parameter: case A and 




50 55 
Snapshot 

Figure 9. Evolution of the circularity parameter for the subsam- 
ple of bound orbits obtained by integrations A (purple curve) and 
B (magenta curve). The integration of the orbits has been per- 
formed on a randomly selected subsample of approximately 100 
galaxies sharing the same infall snapshot, given their initial po- 
sition and velocity. The plot shows the temporal evolution of the 
mean of the circularity distribution for the bound population, as 
a function of snapshot number - from 41 (i.e. 2=1, the infall 
redshift) to 63 (i.e. z = 0). The integration has been split into five 
slots extending roughly 1 Gyr in time; at the beginning of each 
slot (marked by a vertical, dotted line) the position and velocity 
of the satellite was updated to the real value obtained from the 
MS. 
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found in the redshift-zero sample. 



5 DISCUSSION 

In the previous sections we have been drawing conclusions 
from a model which, albeit being state-of-the-art, is well 
known to contain some level of simplification in its treat- 
ment of galaxy formation and evolution; caution is therefore 
needed. 

As discussed in Section 4.4 of iGuo et alj (|201ll ). only for 
galaxies in the stellar mass range 9.5 < log(M*/M0) < 
11 does the model predict u — i colour distribution in 
agreement with observations; at lower (higher) masses the 
galaxies tend to be too red (blue) and the colour distri- 
bution strongly deviates from the reference given by the 
SDSS/DR7 sample. These results were obtained when ap- 
plying the model to both the MS and the Millennium-II 
Simulation l)Bovlan-Kolchin et al.l [20091 ). selecting galaxies 
with log(M*/Af0) > 10 from the former and galaxies with 
log(M*/Mo) ^ 10 from the latter. We have performed 
our analysis on the full galaxy population and after ap- 
plying a mass-cut leaving objects within the above ranges 
(9.5 < log(M*/M ) < 11 and 10 < log(M*/M ) < 11) 
only; the results remain qualitatively the same. 
As already mentioned, galaxies orphans of their dark-matter 
halo, or "type 2s" , were not included in the analysis due to 
the lack of information on their real dynamics. These objects 
are an important ingredient in the model as their presence 
allows the predicted galaxy luminosity function and radial 
number density profile in clusters to match observations; 
in addition to this, they represent a substantial fraction of 
the satellite populations (they account for half of all clus- 
ter members with M* > 1O 1O M0). Including these objects 
in our sample has the immediate effect of increasing the 
median stellar age, while lowering the median specific star 
formation rat^fl; briefly, the sample ages. The effect on the 
anisotropy parameter corresponds to enhanced radial mo- 
tions: from p ~ 0.25 for the full redshift-zero population 
without type 2s, we move to /3 ~ 0.35 when these are in- 
cluded. The relative differences between the blue and red 
population are maintained. The reason for this substantial 
increase in radial anisotropy can be traced back to the ori- 
gin of the orphan population itself; galaxies on more radial 
orbits are more subject to stripping and will lose their dark- 
matter component more easily. It may well be that a dark 
halo is still present in these objects, but that it consists of 
an insufficient number of particles to be regarded as a real 
sub-structure by the halo finders; in these cases, the tra- 
jectory of the most bound particle would still preserve the 
radial feature of the original satellite orbit. 
We have split the original sample of galaxies into two dif- 
ferent populations characterised by different values of the 
u — i colour indicator. As a threshold, we have adopted the 
value where the colour distribution separates into two dis- 
tinct component; for the galaxy sample we selected, this 
occurs at u — i as 2.5 at redshift zero. We have performed 



6 The u—i colour index increases, and correspondingly the sample 
becomes redder, at a given a stellar mass 



the same analysis changing this threshold and making more 
extreme cuts, both at redshift zero and at higher redshifts; 
the results remain coherent with those shown in the previous 
section. We have also used other properties, besides colour, 
to identify the two different galaxy populations, namely spe- 
cific star formation rate and mean stellar age. Not surpris- 
ingly, as all these properties are expected to relate to one 
another, the results have not significantly changed. 
We already mentioned that, according to our results, 
the type of orbit a galaxy is moving on when it last 
enters the cluster environment has major consequences 
on the evolution of the internal properties of the ob- 
ject. The finding must be related to the mechanisms re- 
sponsible for these environment-induced changes affecting 
satellite galaxies moving within a cluster potential. The 
model of GUOll features a more sophisticated and real- 
istic treatment of these effects with respect to earlier at- 
tem pts; indeed, as opposed to previous semi-analytic mod- 
els jBaldrv et al.ll2006l; IWeinmann et al.llioOrl I Wang et"aH 
120071 ; |Pe Lucia fc BlaizotJ 1200/1 ). where the hot gas associ- 
ated to a galaxy is immediately and entirely removed as the 
objects is accreted on to a larger system, in GUOll the 
stripping of gas is performed gradually and modelled to re- 
produce the combined effect of tidal and r am-pressure forces 
(ext ending and integrating the recipes of iFont et all (|2008T ) 
and lWeinmann et al.l (|2010T )). In addition to this, these pro- 
cesses are activated only while the satellite resides within 
the virial radius of the host; this limits the impact of pos- 
sible failures in the FOF algorithm to identify physically 
independent, but spatially close, structures. The improved 
treatment of environmental effects allows GUOll to repro- 
duce, with noticeable accuracy, the radial distribution of 
star-forming, cluster galaxies as found in the SDSS data for 
a large sample of nearby clusters (see their Fig. 3). Albeit its 
successes, the model is still incomplete as it lacks any treat- 
ment for the stripping of the cold-gas component, an effect 
know n to be important in the inner regions of rich clus- 
ters (|Gavazzil 1 19891 ; ISolanes et alj l200ll ; iBoselli fc GavazzH 
120061 ). Even though the overall treatment of environmental 
processes in the model does not entirely encompass the full 
range of mechanisms at work in real clusters, we think our 
results are not significantly affected; if anything, the inclu- 
sion of cold-gas stripping and therefore a more aggressive 
implementation of gas removal could only strengthen the 
differences found in the galaxy properties as a function of 
their initial orbit. 

These speculations are confirmed by the results we found 
when applying a similar analysis to hydrodynamical simu- 
lations, where the co-evolution of dark and baryonic matter 
is followed self-consistently throughout time. We have con- 
sidered few among the most massive and relaxed clusters 
from the HutlQ sample (|Dolag et al.l l20o"9h , a set of high- 
resolution, zoomed cluster simulation with implemented 
cooling, star-formation and feedback processes; we also used 
the results from the cosmological Magneticum Pathfinder 
Simulations^, isolating the most massive clusters in the 
128 Mpc/h-side box ("Box 3" run) as well as in the 896 
Mpc/h-side box ("Box la" run). In all cases, we found 
anisotropy profiles in qualitative agreement with those from 



7 http: / /www. mpa-garching.mpg.de/ ~kdolag/Simulations/ 
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the semi-analytic models in terms of radial behaviour and 
global values for j3. More importantly, we found that simu- 
lated galaxies characterised by either a young stellar popu- 
lation, low u — i colour index or high specific star formation 
rate (> 1 x 1O _11 M0 j/r _1 ) present a systematically lower 
anisotropy parameter in each of the radial bins. 
To summarise, we are confident that the intrinsic limits of 
the semi-analytic model by GUOll do not affect the con- 
clusions drawn in this paper. 

A final comment on our choice regarding the cluster sample: 
we have focussed our analysis on the most massive clus- 
ters found in the MS - objects with virial masses greater 
than 2 x 10 14 Mq - but we emphasise that there is no rea- 
son to expect the results obtained for this sample to ap- 
ply at lower host masses. In fact, studying the properties 
of the dark-m atter subhalo s ident ified in the Millennium- 
II Simulation, iFaltenbacheil <|2010l ) shows that their global 
anisotropy parameter depends both on the host mass and 
on the environment the host is sitting in, with satellites re- 
siding in either massive or isolated clusters moving on more 
markedly radial orbits. On the other hand, IWetzell (|201ll ) 
examines high-resolution, N-body simulation and evaluates 
the orbital parameters of satellites at the time of infall, de- 
fined as the first time the object crosses the virial radius of 
a larger host halo; again, a clear dependence of the initial 
orbit on the host mass is found, with galaxies characterised 
by increasingly radial motions at higher host halo masses. 
We therefore do not expect our results to apply to environ- 
ments other than those of the most massive structures, at 
least not on a quantitative level; also the dependence of the 
internal evolution of the galaxy properties on the initial or- 
bit may not be as strong in smaller objects or groups, where 
ram-pressure stripping is much less efficient a process than 
in rich clusters. 

Besides the results on the link between orbital and inter- 
nal properties of galaxies, we also showed our findings on 
the temporal evolution of the orbital anisotropy itself. We 
found that this quantity increases in time, especially just 
after infall. We note that our definition of infall time may 
have an influence in this analysis, at least on the magnitude 
of the infall anisotropy; half of the substructures hosting 
galaxies become part of the friend-of-friend halo of the clus- 
ter at distances exceeding two virial radii and may therefore 
have an extra tangential component in their initial veloc- 
ity arising from the pull of the matter distribution at the 
outskirts of the host. We found that we can reproduce the 
increase in anisotropy by integrating the orbits of galaxies 
in a spherically-symmetric, smooth mass distribution, when 
the cosmological growth is taken into account. When inves- 
tigating the evolution of the circularity distribution for the 
bound subsample of the integrated orbits, we found very lit- 
tle evolution. This suggests that the relation between the 
anisotropy parameter for a population of galaxies and the 
distribution of their circularities is less naive than a simple 
one-to-one correspondence. 

We register some tension between our results and exist- 
ing obs ervations. Analysing t wo cl uster sets at z w and 
z ~ 0.6. lBiviano fc Poggiantil (|2009T l depict a scenario where 
radially-infalling galaxies progressively turn from being star- 
forming to quiescent, while reducing their anisotropy from 
positive values down to /3 ~ 0. At redshift zero, the re- 
sulting anisotropy profile for the star-forming galaxy sam- 



ple is consistent with more radially biased orbits than those 
characterising the quiescent counterpart. This is equivalent 
to saying that the orbits evolve in parallel with the inter- 
nal properties of the satellites and that the former do not 
significantly impact the latter. Our analysis does, instead, 
suppor t the opposite s c enari o, as thoroughly discussed. Ear - 
lier on, iMahdavi et~aI1 l|l999h and lBiviano fc Katgert] (|2004 ) 
reported of observational findings supporting more radial 
orbits for late-type galaxies than for early-type ones; their 
explanation was, again, in terms of coeval changes in the 
orbital and internal properties of the objects. Other ob- 
servations, as mentioned in Sec. [TJ report of gas-rich, late- 
type spirals in nearby clusters being found o n more tangen- 
tial orbits than ear ly-type, gas-poor spirals (|Dresslerlll986l ; 
ISolanes et al.|[200ll ); this is better reconciled with our find- 
ings, at least those regarding the impact of the satellites' or- 
bits on their internal evolution. We note though, that these 
works agree in finding early-type galaxies (lenticulars and 
ellipticals) on close to "isotropic" orbits. 
Both the observational and the numerical approach come 
with specific limitations and it would be interesting to pin- 
point the origin of the disagreement. We leave this task to 
future works. 
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